Low Temperature Specific Heat of some Quantum Mean Field glassy phases. 
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We investigate analytically the low temperature behavior of the specific heat Cv{T) for a large 
class of quantum disordered models within Mean Field approximation. This includes the vibrational 
modes of a lattice pinned by impurity disorder in the quantum regime, the quantum spherical p- 
spin-glass and a quantum Heisenberg spin glass. We exhibit a general mechanism, common to all 
these models, arising from the so-called marginality condition, responsible for the cancellation of 
the linear and quadratic contributions in T in the specific heat. We thus find for all these models 
the Mean Field result Cv{T) oc T^. 
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I. INTRODUCTION. 

While they have been experimentally observed several 
decades agoi, the anomalous low temperature thermo- 
dynamical properties of disordered and glassy systems 
remain a formidable theoretical issue. In particular, 
measurements of the specific heat Cv{T) in a variety of 
glasses including structural glasses^, disordered crystals^, 
or spin-glasses'^ show a linear behavior Cv{T) oc T at 
low temperature. Such a behavior is often explained by 
the standard two-level systems (TLS) phenomenological 
argument^. Although this TLS argument is very appeal- 
ing, and successfuU in many situations, it appears ex- 
tremely hard to vindicate it from a microscopic descrip- 
tion. 

The computation of the specific heat of a disordered 
system starting from a microscopic hamiltonian is a very 
complicated task. In this respect, important progress 
has been achieved by recent developments in Mean 
Field methods in quantum spin-glasses^iSiLSi and related 
models^, allowing for the description of low lying excita- 
tions in quantum glasses. However, even in this solvable 
limit, the analytical computation of the Cy{T) is still 
intricated and the question whether this TLS argument 
is confirmed or not by Mean Field calculations is still a 
subject of controversjiiSiii. 

In Refii, a quantum extension of the spherical p-spin 
model was studied. In the marginal spin-glass phase, 
characterized by a one step Replica Symmetry Breaking 
(RSB) ansatz together with the marginality condition, 
some indications were found for a linear behavior of the 
specific heat, although its low T behavior was not ex- 
tracted analytically. The authors of Ref.^ have studied a 
Mean Field theory of a SU{2) quantum Heisenberg spin- 
glass. Using a semi-classical expansion in 1/5, with S 
the size of the spins, the specific heat was obtained ana- 
lytically to lowest order in the marginal spin glass phase, 
also described by a one step RSB solution. At this order, 
the linear and quadratic terms of the low T expansion of 
Cy(T) were found to cancel, leading to a cubic behavior, 
Cy{T) cx (these cancellations were found to occur in 
the related model of a quantum Ising spin-glass at the 



lowest order in a similar semi-classical expansion^). The 
expansion to next order appeared to be rather intricated, 
and it was argued^ that the accidental cancellation iden- 
tified to lowest order does not occur to this next order, 
yielding a linear contribution to Cy (T) . A later numerical 
solution of the saddle point equationiS claims instead the 
absence of this linear contribution and a low T behavior 

CyiT) CX r2. 

A class of models for which such Mean Field meth- 
ods have been applied with some success, e.g. to com- 
pute correlation functions2iiS, are disordered elastic sys- 
tems, which cover a wide range of physical situations such 
as charge density wavesi^, electron glasses"'^'*'^^, and flux 
latticesiSiiLiLiS, for which the quantum regime is rele- 
vant. In the elastic limit, where topological defects can be 
neglected, which is for instance the case in the so-called 
Bragg glass phase, these systems have been studied, both 
in the classical^^ and quantum^ limit, using the Gaussian 
variational approximation2ii242& to the replicated Hamil- 
tonian. In this framework, the specific heat has been 
studied both in the classicalSi and quantum regimeili?!?;. 
In the quantum limit, of interest here, the specific heat 
has been computed in a semi classical expansion in pow- 
ers of keeping h/T fixed, similar to the aforemen- 
tioned 1/5 expansion^. At the leading order, the can- 
cellation of the linear and quadratic terms in Cy (T) was 
also obtainediiiSS. But surprisingly, the analysis of the 
next to leading order showed that these cancellations also 
occurS^. In view of these result a^i"'^°i^^ , it is important to 
know whether there is a general property, within this 
Mean Field approach, leading to the cancellation of the 
linear term in Cy(T). 

In this paper, we identify a general mechanism, com- 
mon to all these models, relying on the marginality con- 
dition, which leads to the cancellation of the linear and 
quadratic (in T) contribution to Cy{T) at low T. This 
leads, independently of any semi- classical expansion nor 
numerics, to Cy{T) (X T^. We beficve that this sheds 
light on the Mean Field approximation applied to these 
quantum disordered models. 

The organization of the paper is as follow. In section II, 
we introduce the different models we will be interested 
in, and recall the main properties of the saddle point 
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equations. Section III is devoted to the low T expansion 
itself: we first exhibit the non trivial low T structure of 
the variational equations, therefore extending the previ- 
ous analysis of Ref.^'^'^ at finite T, and then turn to the 
computation of the specific heat. Finally, we draw our 
conclusions in the last section. 



II. MODELS AND MEAN FIELD 
APPROXIMATIONS. 

A. Quantum periodic elastic manifold in a random 
potential (Model I). 

We consider a collection of interacting quantum ob- 
jects of mass M whose position variables are denoted 
Ua{Ri,T). The equilibrium positions Ri in the absence 
of any fluctuations form a perfect lattice of spacing a. 
Interactions result in an elastic tensor <&ct,/3(9) which 
describes the energy associated to small displacements, 
the phonon degrees of freedom. Impurity disorder is 
modeled by a r independent gaussian random poten- 
tial U{x) directly coupled to the local density p{x) = 
5{x — Ri — u{Ri, r)). We will describe systems in the 
weak disorder regime a/Ra <C 1 where Ra is the trans- 
lational correlation length, 6.17. in a Bragg glass phase 
where the condition \ua{Ri,T) — Ua{Ri + a,T)\ <C a holds, 
no dislocation being present. The system at equilibrium 
at temperature T — l/P is described by the partition 
function Z = Tre^l^^ = / DuDUe"'^/'^ with the Hamil- 
tonian H = iJ ph + i?dis : 

1 f nfo)^ 

i^ph=2 J -^+IZ"a(9)$a,/3(<z)M;3(-'7) 



'9 a,l3 

Hdis^ I d'^xU{x)p{x,u{x)) 



(1) 



and its associated Euclidean quantum action in imagi- 
nary time T 

-S[n,u]= dr iU^{q,T)drU^{q,T)-H (2) 

Jo Jq 

where u{q, r) and its conjugated momentum Il{q, r) sat- 
isfy periodic boundary conditions, of period (3h, along 
the r axis. One denotes by = J^^ (^Try integration on 
the first Brillouin zone. We will focus here on the case of 
internal dimension d > 2. For simplicity we illustrate the 
calculation on a isotropic system with ^a,fi{q) = cq'^^ap 
and denote disorder correlations 



The disorder average is performed using the replica trick 
Z'^ — J Due~^"'p^^ and integrating over 11, after some 
manipulation^, one obtains the following replicated ac- 
tion Srcp = Sph + Sdis with: 



^ph = j d'^xdr^Y.i^^Uaf + ^{drUaf 

S'dis = y d'^xdrdr' ^ R{ua{x, r) - Ub{x, t')) 

ab 

R(u) ^ pI^Ak cos{K ■ u) (4) 



K 



Here v = c/M is the pure phonon velocity and Ak = 
/ d'^xe'^^'^ A{x) denote the harmonics of the disorder cor- 
relator at the reciprocal lattice vectors K, and po 
the average areal density. 

Given the complexity of the replicated action 5icp, it 
has been proposed to study it within the Gaussian Varia- 
tional Method (GVM)2ii2iSli. It is implemented by choos- 
ing a trial gaussian action Sq parametrized hy a k x k 
matrix in replica space [q^ujn): 



■^0^9^/ ^'^abil^'^n)u''{q,UJn)u^{-q,~i^n) 
^P'^Jl a,b 

G^^{q,UJn) = Cq^5ab - <yab 



(5) 



U[x) = Q , U{x)U{x') ^ A{x-x') (3) 



which minimizes the variational free energy i^™'' = 
^0 + ^^^'^'^^ ~ So) So, where denotes the free en- 
ergy computed with Sq. In the limit A; ^ 0, we denote 
G{q,ujn) = Gaa{q,ujn) and parametrize Ga^b{q,i^n) by 
G{q,u), where < u < 1, which is ujn independent^. 
Similarly we take Bab{T) = {[ua{x,T) — Ub{x,{))Y) /m 
with B{t) and B{u) which is r independent. The best 
trial Gaussian action (|SJ) is obtained by breaking the 
replica-symmetry (RSB)i2i. A previous analysia^^ re- 
vealed indeed the existence of a breakpoint Uc such that 
ct(m) = a{uc) for u > Uc- In d > 2, where there is a full 
RSB solution, (7{u) is a continuously varying function of 
u for u < Uc- In d = 2, for the single cosine model, there 
is instead a (marginal) one step RSB solution such that 
a{u) = for u < Uc- 

Using the variational approach, it has been shown in 
detailsS^ that the specific h eat i s obtained from the T- 
derivative of internal energy {H) per unit volume, which, 
independently of the RSB scheme, can be written in 
terms of the saddle point solution: 
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cg^ + S + IjoJn) 
cgf2 + E + Mw2 + /(CJ„) 



0h 



dT{F{B{T))- F{B))- / dw{F{B{w)) ~ F{B)) (6) 



r 



where F{X) = V{X) - §V'iX), V{X) = 
-plT,K^Kexp{-XKy2), Wc = f3uc. In ®, the 
quantities entering this expression are determined by 
the variational equations: 



2 r'^'^ 

/K) ^j^j^ dr(l-cosKT))(f' (B(t))-F'(B)) (7) 



1 = -4F"(B)J2(S) 



(8) 



with the definitions 



B=l\- f I 



1 — cos (ijJ„t) 



The breakpoint Wc, in d > 2, is determined by: 



(9) 
(10) 

(11) 



We finally quote the following useful relation, valid a full 
RSB solution, obtained by combining Eq. (jHJ and Eq. 



WcSB + 2J2{^)SJ: = 



(12) 



where d stands for an infinitesimal variation. The equa- 
tion © is the marginality condition, corresponding to 
the vanishing of the replicon eigenvalue, which holds 
automatically in this problem for d > 2. As first no- 
ticed in^ and further investigated in^^-22, the solution of 
the variational equations can be organized in an expan- 
sion in h keeping ph fixed. Expanding any quantity Q 
as Q — J2^=Q^^Qn{0fi), it was shownS that this con- 
dition describes a gapless excitation spectrum caracter- 
ized by the low frequency behavior of the self-energy Q 
Io{uJn) oc \uj„\ + 0(^2) leading to the analytic continua- 
tion /q ('^) ^ ^ where /(w„ —10-1-1-0+) = I'{uj)+il"{u!). 
Thus one expects a power law behavior of the specific 
heat at low T. In the following, we will compute analyt- 
ically the low T expansion of the internal energy l^. 



B. Quantum spherical p-spin glass model (Model 

II). 

We consider a quantum extension of the spherical Jo- 
spin glass model as studied in^, an interacting system of 



N continuous spins s^, 1 < z < A^. This quantum ex- 
tension consists in considering a continuous spin Si as an 
operator associated to a spatial coordinate and introduc- 
ing its conjugated momentum tt^ which satisfies standard 
commutation relations 







'ihSi 



(13) 



The quantum p-spin glass model is then described by the 
following hamiltonian 



_ _ 7r2 ^ 



(14) 



where we denote tt'^ =1t-^, with^ = (tti, ttat) (sim- 
ilarly for s2 and ~s) and impose the spherical constraint 



1 ^ 
iV ^ 



= 1 



(15) 



In (|14|l . the coupling constants t/ii,..,ip s-re random vari- 
ables, independently distributed according to a gaussian 
distribution of zero mean and variance 
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2NP-^ 



(16) 



This model (|14ll is then studied^ using the formalism of 
the quantum action in imaginary time together with 
the use of replicas to implement the average over the 
disorder (|16|l . After some manipulations, one obtains, in 
the limit A^ — > 00 a saddle point equation for the order 
parameter 



Qab{r-r') = -{-^a{T)-^b{T')) 



(17) 



where a, h are replica indices. In the limit fc ^ 0, one de- 
notes Qaa{T) = Idi'T) and parametrizes Qa^b{T) by q{u) 
which is T-indcpendent. Following the authors of Ref.-i, 
we will work with dimensionless quantities, by redefining 
the imaginary time f = Jt /h and Matsubara frequencies 
ibn — hujn/ J (in the following we will drop all hats in 
order to simplify the notations). We also introduce the 
parameter F = jiMJ), which measures the strength 
of quantum fluctuations. The phase diagram of H14I) in 
the F — T plane was foundi to be characterized by a 
line Fc(T) separating a paramagnetic (PM), associated 
to a diagonal matrix Qab(r) = qd{T)Sab, from a spin-glass 
(SG) phase at low T, which we focus on here. The saddle 
point equations describing this SG phase is solved by a 
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one step RSB ansata^-, shown to be exact as in the clas- 
sical case, such that q{u) — for u < m and q{u) = q^A 
for u > m being the breakpoint. The internal energy, 
as a function of the saddle point solution is given bj^: 



'0 



dr{q'Ar)-q^^^) (18) 



with the saddle point equations 



/3 



EK) = I dr(l-cosKr))(9r'(r)-ggrXl9) 





z = -iimq"^^ — 



and the definitions 

^EA 



P^^ + z' + f:{, 



UJ„ 



qd{r) - <ZEA = 



cos LOriT 



n r 

The breakpoint is determined by 



/3 ^ !^ + z' + S(a;„) 



(3m = Xp 



p{xp + 1) 



-p/2 

9ea 



(20) 



(21) 



(22) 



(23) 



Combining (|2U|I and H23|) , one obtains the following useful 
identity 



p{l + Xp) 



2-p 

9ea 



(24) 



As it was noticed in other one step RSB solutions^, one 
obtains a one parameter family of solutions, indexed by 
Xp (or equivalently by the breakpoint m). There are then 
two different ways to determine m. In the statics, m is 
usually determined by minimizing the free energy : this is 
the so-called equilibrium criterion. The excitation spec- 
trum of the equilibrium SG state is gaped, yielding a 
specific heat wich vanishes exponentially at low T. Al- 
ternatively, m is determined by imposing the vanishing 
of the replicon eigenvalue, which leads to the so called 
marginality condition^: 



Xp = p ^ 2 



(25) 



One can show^, using a Keldysh mean field approach and 
performing analytical continuation to imaginary time, 
that the marginality condition H25I) gives indeed the cor- 
rect solution from the dynamical point of view, i.e., if one 
considers, in an infinite system, the large time limit where 
time translational invariance and equilibrium fluctuation 
dissipation theorem hold. Moreover, this marginal value 
of Xp was found to be the only one compatible with 
a gapless excitation spectrum^. In the T = limit. 



it was indeed shown that, in the low frequency limit 
5i(w„) oc \ujn\ + 0{u!^). Therefore, one expects that the 
specific heat of the marginally stable SG state vanishes 
as a power law. In the following, we show how to extract 
analytically the low T behavior of the marginally stable 
SG state. 



C. Quantum SU(N) spin-glass (Model III). 

We consider the Heisenberg quantum spin-glass, de- 
fined by the following hamiltonian^*^ 



H 



1 



NAT 



(26) 



where the original spin symmetry group SU (2) is ex- 
tended to SU {N)S. and the large N limit is taken. These 
J\f spins occupy the sites of a fully connected lattice. In 
(|26|l . the coupling Jy are random variables, indepen- 
dently distributed according a gaussian distribution of 
zero mean and variance 



J' 



(27) 



Using the imaginary time path-integral formalism (we 
will set h—\ from the beginning for this model) , together 
with replicas to implement the average over the disorder 
(|27|l . the model is mapped, in the infinite range limit, 
onto a self-consistent single site problem described by 
the action^: 

SaS = SB~— drdr'Qabir ~ r') S air) ■ S b{r'l28) 

where Sb is the Berry phase imposing the spin commuta- 
tion relations^, together with the self-consistent equation: 



Qabir 



N 



^iSair) ■ S'6(r'))seff 



(29) 



where {..-jsga stands for an average computed with the 
action SeS. Using a bosonic representation of the spin op- 
erator S in terms of Schwinger bosons, Safs = bl^bp — S6af3 
with the constraint ^a^a — SN, this model 1)26(1 
can be described analytically in the limit N —i- oo 
which then constitutes a mean field theory of the fully 
connected model (|26ll where the spins have the sym- 
metry SU{2). In this limit, the original self consis- 
tent equation H29|) reduces to an equation for the boson 
Green's function ^''''(r) = -(T6"(T)6tf'(0)). In the hmit 
fc ^ 0, one parametrizes C/''°(r) by G{t) — g, such that 
lim^^oo G{t) = at T = and, G'^^^'ir) by -g{u), which 
is T-independent. The phase diagram of H26|) in the T — S 
plane has been established in the large N limi^. A line 
Sc{T) separates a paramagnetic phase, described by a 
diagonal matrix in replica space G°'''{t) = SabG{T), and 
where several crossovers were found to occur in the quan- 
tum regime^, from a spin glass phase, which we focus on 
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here. In this SG phase, the saddle point equations are 
solved by a one step RSB ansatz, such that g{u) = for 
u < X and g{u) = g for u > x, x being the breakpoint 
and g — g. The starting point of our computation of the 
specific heat is the expression for the internal energy per 
unit volume^: 



{H) = J^^ driGir) - g)Hg{~T) - gf 

-:^P[x~l)g' 
in terms of the saddle point solution 



S(ii/„) = r / dr(e-*'^"^ - 1) 
X ({g{T)~gf{g{-T)-g) + g^ 

-i- 

Jg^ \Q 



(30) 



(31) 



(32) 



with the definitions 



g = s + g{T = o-) 



g{ivn) = 



1 



- f - E(zz.„) 



(33) 



where Vn is a bosonic Matsubara frequency. Similarly to 
the spherical p-spin model H19I23|I . one obtains a one pa- 
rameter family of solutions, parametrized by 0, or equiv- 
alently by the breakpoint x. Here also, if one chooses the 
equilibrium criterion, the excitation spectrum is gapped. 
Instead, if one imposes the vanishing of the replicon 
eigenvalue, one obtains the marginality condition^: 



e = e 



1 



R 



V3 



(34) 



Using an expansion in 1/S - similar to the semi-classical 
expansion for the elastic manifold^ -, it has been shown 
explicitly^ that the marginality condition (|34|l is the only 
one compatible with a gapless excitation spectrum, such 
that t"{uj) oc uj + O(w^), where E(i/„ -iu + 0+) = 
S(w) -|- iE"(w). Although the connexion between this 
Matsubara formalism using the marginality condition 
(|34|l and the true hamiltonian dynamics in real time at 
the Mean Field level has not yet been established for the 
present case, the study of similar one step RSB solutions 
for which this connexion has been don o^i^^ , suggests that 
(|34|l gives indeed the correct solution from the dynamical 
point of view. On physical grounds, in the present case 
of the ordered phase of a quantum SG with continuous 
symmetry, this choice seems also natural as one indeed 
expects a gapless excitation spectrum^. 

Although the specific heat in the equilibrium SG state 
vanishes exponentially at low T, we focus here on the low 
T behavior of (T) in the marginal state H34|) , which is 
instead expected to vanish as a power law^. 



III. LOW TEMPERATURE ANALYSIS. 

In this section, we compute the low temperature ex- 
pansion of the internal energy for the different models 
presented before, from which we directly obtain the spe- 
cific heat Cv{T). To do so, we start by deriving some 
general identities, which form the background of our anal- 
ysis. The finite temperature behavior of the internal en- 
ergy requires the low temperature expansion of the sad- 
dle point equations, which we then present in details. For 
that purpose, we will use the notation, for any quantity 
of interest Q,Q = Y.n=o T"""?^"^- We finally turn to the 
behavior of Cy{T) in the last paragraph of this section. 



A. General properties. 

We first focus on the two first problems evocated here 
(|1I14|) . which show, formally, a strong similarity. In par- 
ticular, at variance with Model III H2()l) . these two sys- 
tems exhibit a two point Green's function (|10l 1221) wich 
is invariant under the transformation a;„ —^ —ujn- We 
focus on the low T behavior of integrals over imaginary 
time r which enter both the variational equation and the 
computation of the internal energy. We will use the no- 
tations of the disordered elastic hamiltonian H7I11|I . the 
transposition to the p-spin model being straightforward. 

We first suppose that E and I{lu„), as a function of w„ 
are independent of T, i.e. E^") = 0, Vn > and similarly 
for I{ujn) with /(w„) ~ |a;„| -f 0(0;^). Then, as shown in 
Appendix one has the low temperature expansion, for 
any function 7i(X) that can be expanded as a power serie 
around 0, n{X) = J2Zo ^kX^ ■ 



dr cos {LUnT){H{B{T))-H{B)) ^C'*+0{T'^) (35) 



where here, and in the following, C"** stands for a 
(generic) quantity independent of T (it may however de- 
pend on the Matsubara frequency). Indeed, the gapless 
structure of the spectral function suggests that only even 
powers of T should enter this expansion H35|) . But as 
shown in details in the Appendix ^ one can check ex- 
plicitly that the term cx vanishes. This property H35() 
is rather independent of the structure of the variational 
equations, only requiring I" {to) oc w at low frequency. 

For the other model discussed here (|2SJ), one could not 
show a so general statement (I35II involving any func- 
tion H, due to the absence of the symmetry if„ 
—ivn of the Green's function g{iVn) H33|l , which ren- 
ders the calculations more subtle. Nevertheless, with 
the same assumptions as above that the complete self- 
energy Jg/9 + E(w„) is independent of T, with 
Yj{ivn) oc \vn\ + 0{vf^), for the particular integrals in- 
volved here, one can show (see Appendix^ for more de- 
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tails): 
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driiGir) - gf{g{-T) - gf - g^) = C^' + 0{T^) (36) 



/3 



dre— "-((^(r) - gf{Q{~T) ~ g) + g') = C^' + 0{T^) 



These properties are very useful tools to investigate the 
low T behaviors of physical quantities in these models. 



B. Variational equations. 

We now turn to the low T expansion of the saddle point 
equations for the three different models. 



1. Model I. 

Most of the properties presented here, and their ex- 
tension to the spin-glass models (|14l I26|l , have been sug- 
gested by the expansion in powers of h, keeping [3h fixed. 
We will shortly remind here the main features of the semi 
classical expansion of the variational equations {TUHJ. We 
use the notation, for any quantity Q = X^fc^o ^^Qk{PK). 
Although at lowest order, the complete solution of the 
variational equations^ shows that Sq ^^nd Io{ijJn) are in- 
dependent of T, Si and Ii{ujn) become T dependent at 
the next orde»2^. At low temperature the following struc- 
ture was explicitly obtained 



= +{^\^^^ +0{{T/Kf) 



(37) 



T 



/f^(cc>„) = -(l-<5„,o)Si- 



(2) 



(38) 



such that the finite temperature corrections of Si -I- 
Ii{oJn) are oc for a;„ ^ 0. It was also noticed that 
the peculiar term H38I) . once inserted in B 0, generates 
odd powers of T, the lowest one being T^, in the low T 
expansion at higher order in the expansion in h. We ex- 
tend here these properties l|37l I38|l independently on the 
semi-classical approximation. We show indeed that the 
following expansion, up to order 0{T^): 

S = S(") + T^S(2) + T^S(^^ (39) 
/(^„) = /W(^„) - (1 - (5„,o)(t2s(2) +t3s(3)) 

with /(°^(cj„) (X \ujn\ + O(w^), is a consistent solution 
of the variational equations (|7I8|I . To do so, we compute 
the low T expansion of the r.h.s of the equation for I{uJn) 
given the forms H39|l . We thus need an extension of the 
general property H35|l . when S and I{uJn) have the form 
given in Eq. (|39|l . As shown in Appendix ^ , the low 



T behavior of such integrals over r (I35() are in that case 
given by 



/3h 



dT cos {uJnT)iniB{T)) ~n{B)) = C 



(40) 



+2M„,o(r2s(2) -f T3s(3))J2(s(*'))H'(S(")) +C'(T'') 

If one applies this general formula H4l)|) with H ~ 2V'{B) 
to the r.h.s of Eq. one obtains up to order 0{T'^): 

-J dT{l~cos{oj^T)){V'{B{T))-V'{B))=C^' (41) 

+4(1- j„,o)(r's(2) +t3s(3))J2(s("))J2(s(o))f"(b(o)) 

^ - (1 - (5„^o)(T2s(2) + T^sf^)) (42) 

where, in the last line, we have used the marginality con- 
dition (|SJ| at T = 0. This relation (|41|l thus shows explic- 
itly the consistency of the low T expansion (|39|) proposed 
for the exact solution of the variational equations. Impor- 
tantly, although the general property shown above H40() 
holds independently of the saddle point equation, a so- 
lution such as H39|l is consistent provided the marginality 
condition holds. 



2. Model II. 

Inspired by the previous analysis, we show that a con- 
sistent solution of the variational equations for the quan- 
tum p-spin model (|19l I2U|) is given, up to order 0{T'*') 
by: 

z' = z'(")+T2z'(2)+rV(3) (43) 
S(c.„) = S(«)(a;„) - (1 - <5„,o)(T^^'(2) + T^z'^'^) 

with S^°)(aj„) c>c \ujn\ + 0{tJ^). This is shown by using 
the general low T expansion, an extension of H40() to the 
present case (|43|) : 



Jo ' 



dr cos {iUnr){H{qd{T)) - H{qEA)) = 



^S„.oiT'z'^'^ + T'z"^'^)^^ + 0{T^) (44) 

Thus applying (gH) with n{X) = {p/2)XP-'^ yields the 
low T expansion of the r.h.s of the variational equation 
(P|l up to order ©(T^) 
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2 Jo 

p{p- 1 



dr(l~cos(c.„r))(gr'M-9gA) = C^' 



■(l-<5„o)(T^^'(2)+rV(3))i|A^ 



^ , (i-^,,o)(rV(2) + r3^'(3)) (45) 

I + Xp 

where we have used, in the last line, the identity H24|l . 
Thus, one sees on (|45ll that the expression given in Eq. 
(|43|l is a consistent solution of (|19|l provided Xp = p — 2, 
i.e. the solution is marginally stable H25|) . 
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3. Model III. 

For this model, we show that, similarly to the two pre- 
vious ones, a consistent solution of the variational equa- 
tions is given, up to order 0{T^) by 



(46) 



J 



t{^v^) = t^'\iv^) + |(<5„,o - l)(T^g(2) +t3^(3)) 

with oc \i>n\ + 0{v'^). To show the consistency 

of this solution (|46(l . we perform the low temperature 
expansion of the r.h.s of the equation for E(i^'„) Ij^lf) 
given H46(l . One obtains up to order 0{T^): 

J2 f dr(e— - l)[{g{T) - g)^(g(-r) - .g) + 
Jo 

= C"' -I- {5nfi - l)3ej (r2g(2) 4_ T3g(3)) (47) 

Again, this expansion (|47|l shows that the structure ex- 
hibited in H46|) is a consistent solution of the variational 
equation provided 36 = 1/9, i.e. 9 = 9_r = l/\/3, 
the marginality condition (|34fl . 



C. Specific heat : low temperature expansion. 

We now turn to the computation of the low tempera- 
ture behavior of the specific heat. 

1. Model I. 

Our starting point is the expression for the variational 
internal energy (H) given in jnj. We first analyse the low 
temperature behavior of the first term in ^ , namely the 
sum over Matsubara frequencies. Importantly, we notice 
that, in this sum, the contribution of the mode uJn — is 
independent of the peculiar structure of /(w„) exhibited 
in H38|l : this allows to avoid ambiguities coming from 
the analytic continuation of such a term cx (1 — Sn.o) in 
(|38|l . We can thus safely transform this discrete sum in 
an integral: 

f cg2 + E + /(L^„) 



cg2 + S + Mu;l + /(w„) 



(48) 



+00 



duj 



(49) 



PDOs(t^) 



I, {cq^ _ -i-c^2 + s + I'iuj)^ + {I"{uj)Y) 

where Pdos('^) is the density of states. Using that 
Puos{^) = Pdos(w) -I- T'^p'^lsiuj) (123, together with the 
low frequency behavior pixls oc w^, one obtains straigh- 
forwardly 

cg2 + S + /(w„) 



-y 



P^Jg C(?2 + E + Afc^,2 + /(c^„) 



C-t+0{T^) (50) 



One performs the low T expansion of the second term 
in (O, the integral over r, using the property H40|l which 
holds for the solution of the variational equations we have 
found 



I3h 



dT{F[B{T)) - F{B)) = 



(51) 



One finally tackles the analysis of the last term in ((HJ, 
the integral over w, by using the property that B{w) is 
independent of T for w < w^. 

/ dw{F{B{w)) - F{B)) ^ C"* 
Jo 

+wi°\T^B^^^ + r3s(3))^'(5(o)) + 0(j^4) (52) 

This property (|52|l is then obvious for a full RSB solution 
where B{w~) = B and was shown^S to hold also for the 
marginal one step RSB solution in d = 2. 

Using the consequence of the marginality condition 
(|12|l . one thus sees clearly that the quadratic and cu- 
bic contributions to this low T behavior cancel between 
Eq. (inUl and Eq. And this yields 



{H) = c"' -i-e'(r^) 



(53) 



This results in the low temperature behavior of the spe- 
cific heat: 



a(T) cxT3 + C'(t4) 



(54) 



Although the coefficient of this cubic term is very hard 
to extract by the method presented here, it has been 
explictly computed at the lowest order in the aforemen- 
tioned semi-classical expansion and it was found to be 



non zero 



2. Model II. 



We analyse the low T behavior of the internal energy 
within the variational method H18(l of the quantum p-spin 
model. We treat separately the two lines of H18f) . Using 
(|44|l . one obtains 



(55) 



Using this expansion and the identity (|24l) one obtains 
that the quadratic and cubic terms in the low T expan- 
sion of the first line in (|18(l cancel provided Xp = p — 2 

+ dr{qT\r) - gg^') = C^' + 0[T^) (56) 
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We now focus on the second line of (|18() . Using (|44|l , one 
has immediately 



IV. CONCLUSION. 



p + 2 



Jo 



(57) 



Finally combining (|23|l with (|24() , one obtains the expan- 
sion of the remaining term in (|18|l up to order 0{T^): 



p+2 p _ p{p + 2) Z^(2)r2 + Z^(3)y3 



(58) 

Collecting the different contributions to (H) H551 15715^ . 
one obtains 



(59) 



which leads to the low temperature behavior of the spe- 
cific heat 



3. Model III. 



(60) 



The low T expansion of the internal energy of the 
Heisenberg spin-glass model (|3U|) is performed using the 
relation derived in Appendix IbI which . given the solution 
of the variational equations we have found, yield 



(61) 



The expansion of the last term in (|30|l is straighforwardly 
computed using the relation (|32|l 

--^5*°^ - (5(2'r2 + 5(3) j^3) _^ 0(y4) (g2) 

Thus combining Eq. (|61|l and Eq. H62() . one sees that 
the quadratic and cubic terms come with a prefactor 
30 — I/O, which thus vanishes for the marginally sta- 
ble solution, corresponding to 9 = 6^ = 1/V3. This 
yields the low temperature behavior of the specific heat 



(63) 



For this model too, the amplitude of the cubic term has 
been computed in a 1/S, semi-classical, expansion, and 
found to be non-zero^. 



To sum up, we have computed the low temperature 
specific heat of a rather wide class of quantum disor- 
dered systems with continuous degrees of freedom, using 
a Mean Field approximation, including disordered elas- 
tic systems in d > 2, the spherical p-spin-glass and the 
Heisenberg spin-glass. For all these models, we have ob- 
tained that the Mean Field approximation yields the low 
T behavior C„(T) cx If ^WM^ . For the Heisenberg 
spin-glass model, the cancellation of the linear term in 
Cy{T) obtained here is in agreement with the numeri- 
cal solution of the saddle point equation of Rcf.^*'. And 
the non trivial structure of the saddle point solution elu- 
cidated here (|46|l could help to clarify numerically the 
status of the quadratic contribution to Cy{T) obtained 
in Ref.ia. 

This Mean Field result Cy (T) cx is at variance with 
the linear behavior commonly expected from two-level 
systems^. As we have shown, the cancellation of the lin- 
ear and quadratic contributions to Cy{T) strongly relies 
upon the marginality condition. And it is worthwhile 
to notice that the physical picture associated to this 
marginal stability criterion, which enforces the (quan- 
tum) dynamics along the flat directions of the free en- 
ergy landscape, seems qualitatively different from the ar- 
gument steming from TLS*. 

For the case of manifolds, our result, within the Gaus- 
sian Variational Approximation also applies to non pe- 
riodic elastic structures, e.g. domain walls, when they 
can be solved by a full RSB ansatz (or its limiting case 
of a marginal one step RSB). And although this Mean 
Field approach is always an approximation for the pe- 
riodic case, it becomes exact for the non periodic one, 
in the limit where the number of components of the dis- 
placement field becomes infinite, as it is for the spherical 
p-spin-glass model (|14ll . in the limit iV — > cx), or for the 
Heisenberg spin-glass model (|26|l when both N,Af ^ 00. 
An outstanding question remains to know whether and 
how this result is modified away from mean field, which 
clearly deserves further numerical and analytical investi- 
gations. 
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1. A first stage with multiple Matsubara sums. 

To begin with, we restrict ourselves to the case where S 
and liuin) as a function of ujn do not depend on temper- 
ature T, for all LOn including the mode w„ = 0. We also 
assume the low frequency behavior I(oJn) ~ I'^nl + C'(w^) 
(UI). Here we are interested in the low temperature ex- 
pansion of the following quantity H35|) : 



J dr {n{B{T)) - n{B) 



(Al) 



As we will see, the first non vanishing finite tempera- 
ture correction is a priori of order : we show that 
this contribution in fact cancels for any (smooth enough) 
function 7i(A"). To show this cancellation, we show the 
following property, for any integer m: 



(3h 



We introduce de notation 

K{lj^) = MMiul + S + /(lj„)) (A3) 

Inserting the definitions of S Q and B{t) ((TU|1 in ||A2|1 . 

one obtains 



I3h 



(i?(T)" - B") = (A4) 



dr 



{2hY 

^'^ ' ni,..,n,„fc=l 

xif(cj„j..is:(w„jii:(w„,^j..i^(w„„) 



H Y^^~^f^^m \ drcos(u;„,T)..cos(cj„,r) 

„ "'0 



where = m\/{{m — k)\kl). Performing the integral 
over r, using the property of parity K(LUn) — ii'(— w„), 
this expression can be written as 



13 h 



dr 



[B{n 



-2fiJi(E)mB'"-^ (A5) 



xK{u;ni + ...+u;nk-i)B 



ni ,..,nfc_i 



APPENDIX A: LOW TEMPERATURE 
EXPANSION FOR MODEL I AND MODEL II 
DETAILED CALCULATIONS. 



In this appendix, we show some general properties of 
the low temperature expansion of multiple sums over 
Matsubara frequencies. We will use here the notations 
of the elastic manifold (|7I11|I (the extension to the quan- 
tum p-spin model being straighforward). 



We now focus on the low temperature expansion of the 
multiple sum over Matsubara frequencies. In that pur- 
pose, we use the spectral representation of the Green's 
funtion : 



1 



cg2 + Mw2 + E + J(w„) 



Aiq,io) 



— A{q,L0)- 

TT lUJr, — UJ 



(A6) 
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where the spectral function A{q, u) is the imaginary part 
of the retarded function and we remind —iuj + 

0+) = -f'(t^) + il"{uj)- All the temperature dependence 
is then contained in the different Bose factors: 



1 



{nf 
1 



r E 

ni ,..,Tifc_i 



dei..dekA{€i)..A{ek) 
1 



(A7) 



E 



ni,..,..„-i 

^ J dei..dekA{ei)..A{ek){fB{ei)- fsiek)) (A8) 

^2)).. 

(A9) 



x(/B(e2) - /s(efe - ei))(/B(e3) - /s(efe - d 
^, {fsiek-i) - /s(efe - ei - £2 - .. - e/c-2)) 



where .4(ti') = v4(<;, w) and /s(e) the Bose factor. This 
expression ()A7|) has a very interesting structure which 
allows us to extract simply the term of order T^. In- 
deed, considering the low temperature expansion of the 
following term (which is the analogous of a Sommcrfeld 
expansion in the fermionic case) for any function Ti.{x) 
whith n{fi) = 

/ deH{e)fB{e-^l) = - T deH{e) 

J J ~oo 

iy^H\f,) + 0{T^) (AlO) 

notice of course that the asumption 7i (/i = 0) is of course 
crucial here IjAlOp . This expansion (jAlOp allows us to 
write formally the Bose factor, when inserted in an inte- 
gral over frequency e: 

(All) 

where 9{x) is the step function {9{x) — 1 for a; > 0, 
9{x) = if a; < 0), and the notation stands for a 
derivative acting on the function which enters multiplica- 
tively with fsie) the integral over e. This form ljAll|) is 
very suitable to extract the coefficient of the term of or- 
der 0{T^) in (|A7|I . Indeed, using HA11() . one discovers 
in (|A7|) , that only the terms where the Bose factors have 
only one frequency in their argument do contribute to or- 
der 0{T^). If we expand, for instance /^(efc — ei — €2) in 
(|A7I) one obtains, for fc > 3 



1 



1 



dei..defc(5(-efe + ei + e2)de^ A{ei)..A{ek) 



{(3hY 3 7r'= 

X (-0(-ei) + 6{-ek)){-e{-e2) + e{-eu + ei)) 
{-e{-ek-i) + 0{-ek + ei + £2 + .. + efc-2)) 



Cfc - (ei + .. + Cfc-i) 

Because of the function S{ek — ei 
the derivative of the terms d^^{—9{ 



(A12) 



-62), one must take 
-e2)+f(-efc+ei) = 



—^{—i-k + ei) (which is a priori the only non vanishing 
term), which yields 



1 1 



dei..defe(5(-efe + ei + e2)'^(-efc + ei) 

(A13) 



[fihY 3 -K^ 

x^(ei)..^(efe)(-^(-ei) + ^(-efc))- 
{-d{-tk-i) + ^(-efc + ei + £2 + .. + efc-2)) 



e/c - (ei + .. + e/c-i) 



= 



We can easily generalise this mechanism to Bose fac- 
tors having more than 3 frequencies in their arguments. 
The case of Bose factors having 2 frequencies in their 
argument has to be treated separately. This one, i.e. 
/B(efc - ei) m yields 

1 TT^ 1 /" r 
JpEj^Y^ J ^'^^■■'^'^kS{-ek + ei)d^^ Aiei)..A{ek) 

X i-d{-e,) + ei-ekm-9{-e^) + 0(-efe + 61 + e^).. 
{-0{-ek-i) + 0{-ek + ei + £2 + .. + £^-2))] ^^^^^ 



fife - (ei + .. + Efe-i) 



Here again, one has to take the derivative de^ {—9{—ei) + 
9{-tk)) = -(5(-efc), which yields 

1 TT^ 1 f 

2~^~k / dei..dek5{-ek + ei)S{-ek)A{ei)..A{ek) 



{fihY 3 TT 
x(-^?(-e3) + ^(-efc + ei+e2)).. 
{-9{-tk-i) + 9{-ek + ei + £2 + .. + efe-2)) 



Cfc - (ei + .. + ek-i) 



= 



(A15) 



where we have used ^1(0) = 0. We conclude from this 
analysis that only the Bose factors having only one fre- 
quency in their argument, such as /s(ei), /_B(e2).. do con- 
tribute to the coefficient of order in IIA7p . Notice that 
they all give the same contribution, thus the low temper- 
ature expansion of ljA7p can be written for /c > 3 

TFT E K{u^,)..K{LO^^_,)K{LOn,+ ■■^i^n,^.,) 



ni,..,nfc_i 



^ :^^'(0)Jfc-l+O(r') 



^^^^ j dei..dekA{ei)..A{ek){-9{-ei)+9{-ek)) 

X (-^^(-£2) + 9{-ek + ei))(-0(-e.3) + 9{-€u + ei + ea)) 
{-9{-ek-i) + 9{-ek + 61+62 + .. + 6k-2)) 



e/c - (ei + .. + Cfe-i) 



(A16) 



where we have used A{~6) — ~A{6) to treat the term 
/^(efc) in ljA7p . For the particular case fc = 2, one has 



2 A'iO) 



(3n 



J2 + j^m)lA'{Q) + 0{l3-^) 



de^+0{T^) 



(A17) 
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where we have used / deA{e)/e = 7rJi(S). Finally, we 
need the low T expansion oi B 



Using these behaviors, one obtains: 



B = + 2?i f ^^'(0) + 0{T*) (A18) 
I 



I3h 



dr 



(2ft)Vi(S])yl'(0) Jm(m - 1)(B(")) 



in — l 



(A19) 



k=3 



k=3 



First, we notice that the two first terms cx just cancel isolating the term cx (1 — Sn,o) in (|39|) : 
and moreover 

E + /K) = E + C(l-(5„,o)+/K) (A22) 

2_^(2n) (-1) c^[k-A [U)jk-i[B'^ ') E = +r2E(2) +t^e(^) +o(r**) 

m-l C = -T2e(2) _ 2^3^.(3) _^ Q(^rj.4^ 

+ J2 {2h)H-l)'Ci{m - k){2h)^A' iO) Jk{B(°^r~'^-') />„) ^ /0(^^) ^ 0(j.4) 

fc=3 

where I{ujn) is defined such that there is no more am- 
biguity concerning its analytical continuation. We want 
to study the T dependence of integrals over imaginary 
time such as (IA5|) when the solution of the variational 
equations are of the form H39|) . In order to extract the 
quadratic and cubic terms in this expansion, we follow 
the previous analysis, except that, here, the mode a;„ = 
, . has to be treated separately H39|l . Notice however that 
this ambiguity does not exist for the computation of -B(t) 



: ^^'(0)( - {2hfiCl,J2{B^''^r-^ 

■m 

-Y,{-l)H2KfJk-i{B^''^r-\kCt 



k=4 



-(m-fc + l)C^i)^ 

j rn(m — l)(m — 2) tt 



^ 3 



This identity HA20p combined with (|A19|) yields finally 
to the announced property 



(fTn|l . We analyse the low T behavior of (|A5|) in the fol- 
lowing way: 



dr 



(b(7 



(3h 



dr 



{Birf 



B" 



fi Jo 



dr 



B" 



-0{T-^) (A21) 



B" 



as announced in the text (|35|l . 



+l3h{B''' 

S=|^Ji(S + C + MloI + /(c^„)) 



(A23) 



(A24) 



2. Handling the peculiar term oc (1 — 5„,o)- 

We now consider the extension of this property (|A21|I 
to the case where E -t- /(w„) now depend on T, and are 
of the form We want to follow the same steps as 

previously ljA2ip . and use the spectral representation of 
the Green function. Therefore, we need in principle to 
know the analytical continuation of the term oc (1 — 6nfl) 
in H39|) . In order to avoid this ambiguity, we start by 



Using the previous property (|A21(I . one has simply 

-J dr (b(t)™ - j = C'* + 0{T^) (A25) 

and the leading corrections at finite temperature are then 
given by the second term in (|A23|) : 



(3{E 



(A26) 



C 



B' 



i — k 
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since 5^°^ = _B(°\ and more generally, one can write, 
for any tOn and any function Ti.{X) the generalization of 
ijMTIl to this case ^ 



I3h 



dr cos 



Kr) (H{B{T))-n{B) 



(A27) 



as quoted in the text H4U|) . 



APPENDIX B: HEISENBERG SPIN GLASS 
LOW TEMPERATURE EXPANSION. 



in the previous models (|A7|I : 



1 /■ ^ 

i=i 



(/i3(e2) - /i3(e4 - ei)) 



(/s(e3) - + £2 - £4)) 



£4 + £3 + £2 - £1 



(B5) 



Under this form, we analyse straighforwardly the low 
temperature behavior, using the property demonstrated 
previously that only Bose factors with one frequency in 
their argument do contribute to this sum. Using the ex- 
pansion ljA10|) . we obtain up to order ©(T^): 



1. General properties. 

We present here the detailed analysis of the low tem- 
perature behavior of integrals over imaginary time of the 
form (|Sn|) : 



driiGir) - gnOi-T) ~ ^ g'] (Bl) 



We first consider the case where the Green's function 
Giwn), as a function of zf„ is independent of T. Similarly 
to the analysis performed in Appendix^ we develop the 
integrand in (|B1|I and perform the integral over r : we 
are then left we multiple sums over Matsubara frequen- 
cies. Depending on the number of Green's functions they 
involve, the integral over r generates 4 different types of 
terms: 



1=^2"^ (B2) 



iy\ ,tl/2 ■,'11^3 



l2 = -2g-^ ^ [g{iiyi)G{ii^2)G{ii'i+ii'2) 
+G{ivi)g{w2)Q{ivi - iv2)\ 



■ivi) 



/3 

Ti = -Ag^g{iiyn - 0) 

We follow the standard analysis and use the spectral rep- 
resentation of the Green's function to handle these terms: 

/OO 7 -1 
—p{u;)-. (B3) 
-00 ll^n - UJ 

= -Im^(w„^tj-f iO+) (B4) 



Performig then the sums over the Matsubara frequen- 
cies, we are left with the same kind of structure as found 



1 r ^ 

^1 = delde2de3de4l[p{e^){-e{-el) + 9{-ei)) 



(-^(-€3) + g(-£l-£2 + £4)) 
£4 + £3 + £2 - £1 



+T^^p'(0) J deide2d€3p{ei)p{€2)p{€3) 

xM(-.,) + ^(-£3))^^fc^Mz£l±M (B6) 

£3 + £2 - £1 

where we have used the relation (which is valid although 
p(e) is not an odd function): 



deide2d£3p(£i)p(£2)(p(£3) + p(-£3))(^'(-£i) - ^'(-£3)) 
M(-£2)+0(-£l+£3)) 



£3 + £2 - £1 







(B7) 



Indeed, one shows this relation (|B7|1 by noticing that 
the step functions reduce the interval of integration to 
£3 > 0, £2 > 0, ei < and £3 < 0, £2 < 0, £1 > 0. A simple 
permutation of the integration variables then lead to the 
relation ||B7)| . 

The analysis of X2 requires the low T expansion of g: 



= g'^°^ +T^^p'{0) + OiT'') 



(B8) 



X2 can be written as 



I2 = 2g^ J deide2de3p{ei)p{e2)pie3) 

„r (/B(£i) - /b(£3)(/b(£2) - /b(£3 ~ £i))) 

^ 1 

£1 + £2 - £3 

, (/i3(£l) - /i3(£3)(/s(£2) - /b(£i - £3)))^ 



(B9) 



£1 - £2 - £3 

from which we obtain the low T expansion up to order 
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I2 = C^' - T^^p'iO) I deide2de3p{ei)p{e2)p{e3) 
( -9(^62) + ei-e^ + 63)) 



, 25(0) p'{0) 



37r /32 



deide2p(ei)p(e2)[4 



£3 + £2 - ei 

ei)-0(-e2) 



£2 - ei 



(BIO) 



ei + £2 

One obtains in a similar way the expansion of I3 , writing 



^3 = ^1 deide2p{€i)p{€2)[4: 



+2 



/B(e2) - /i3(ei) 



ei - £2 



/^(ei) - /B(-e2)i 



ei + £2 

form which we obtain the low T expansion: 



(Bll) 



I3 = C''+T^^-p'{0) J de,de2p{ei)p{e2lBU) 



£2 - £1 £1 + £2 

+THn[g(°^Yp'{Q)g{iVn = 0) 



(B13) 



where we have used / dep{t)/e = ■nQ{iVn = 0). As Q{ii'n) 
does not depend on T, one obtains the expansion of X4 
up to order: 

X4 = C^* - T24^[5(o)]2p'(0)^(z:/„ = 0) + 0(T4) (B14) 

Finally, collecting the quadratic contributions in ljB6l 
iBlOl iBlll IbT4| yields 



2= / dT[{g{T)-gf{g{-T)-gf-g']^C^' + 0{T^) 
Jo 

(B15) 

as announced in the text (|36|l . 



Using the same method exposed previously, we study 
the T- dependence of the integral over r which enters the 
variational equation of this model 1)3 l|l : 



J= I dre— "-[(5(r)-.g)^(t;(-r)-.9) + g'^](B16) 

'0 

As previously, after we have performed the integrals over 
T in l|B16|l . we have to handle exactly the same inte- 
grals as in l2,^3,^i- The mechanism of cancellation of 
the quadratic term is again at work here (notice however 
that, given that the integrand is here a polynom of de- 
gree 3, the terms like (|B7|) are not present here). This 
yields 



(B17) 



as given in the text H3()|l . 

2. Handling the peculiar term oc (1 — Snfi)- 

We generalize these properties to the case where the 
solution of the variational equations is of the form shown 
in Eq. (|46|l . We use here the same strategy as presented 
in the Appendix ^ We isolate the mode w„ = as 
in l|J23l IJ24|) . and given that (^(r) - g) or (^(-r) - 
g) do not depend on this ivn = mode, one obtains 
straightforwardly: 

I = / dr[{g{T) - gf{g{~T) - gf - = C^* 

JO 



(B18) 



3= / dTe-^^-^mT)-gY{g{-T)-g) + g'] = C^' 
Jo 

as anounced in the text (|47I61|I . 



